Repulsive Fermi gas in a harmonic trap: Ferromagnetism and spin textures 
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We study ferromagnetism in a repulsively interacting two-component Fermi gas in a harmonic 
trap. Within a local density approximation, the two components phase-separate beyond a critical 
interaction strength, with one species having a higher density at the trap center. We discuss several 
easily observable experimental signatures of this transition. The mean field release energy, its 
separate kinetic and interaction contributions, as well as the potential energy, all depend on the 
interaction strength and contain a sharp signature of this transition. In addition, the conversion 
rate of atoms to molecules, arising from three-body collisions, peaks at an interaction strength just 
beyond the ferromagnetic transition point. We then go beyond the local density approximation, and 
derive an energy functional which includes a term that depends on the local magnetization gradient 
and acts as a 'surface tension'. Using this energy functional, we numerically study the energetics of 
some candidate spin textures which may be stabilized in a harmonic trapping potential at zero net 
magnetization. We find that a hedgehog state has a lower energy than an 'in-out' domain wall state 
in an isotropic trap. Upon inclusion of trap anisotropy we find that the hedgehog magnetization 
profile gets distorted due to the surface tension term, this distortion being more apparent for small 
atom numbers. We estimate that the magnetic dipole interaction does not play a significant role in 
this system. We consider possible implications for experiments on trapped ^Li and '"'K gases. 
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I. INTRODUCTION 

In recent years, a series of beautiful experiments have 
shown that a gas of two-component fermions interacting 
via a Feshbach resonance exhibits a superfluid state at 
low temperature [1~4]. An exciting new direction for cold 
atoms experiments is the study of ferromagnetism aris- 
ing from repulsive interactions in a two-component Fermi 
gas. Such a ferromagnetic 'Stoner instability' [5] occurs, 
within mean field theory of a homogeneous Fermi gas at 
zero temperature, when the (repulsive) s-wave scatter- 
ing length, 05, between two spin states is large enough 
that kpo-s > i'/2, where hkp is the Fermi momentum of 
the gas. This condition can be satisfied upon tuning as 
to large positive values near a Feshbach resonance pro- 
vided the system stays stable for a sufficiently long time. 
Since ferromagnetism arises from two-body interactions, 
whereas atom loss due to Feshbach molecule formation is 
because of three-body collisions [6] , there may be a range 
of densities where the lifetime is long enough to reach the 
ferromagnetic state. 

The suggestion that ferromagnetism may be achieved, 
as a metastable state, in cold Fermi gases is not new. 
Salasnich and co-workers [7] have studied the mean field 
theory of a harmonically trapped Fermi gas with repul- 
sive interactions and found that this should lead to phase 
separation between the two spin species if the net magne- 
tization is zero. A similar study was carried out by Sogo 
and Yabu [8] allowing for non-zero spontaneous magne- 
tization. Duine and MacDonald [9] later showed that 
the ferromagnetic transition in a homogeneous Fermi gas 
changes from a continuous to a first order transition at 
low enough temperatures upon going beyond mean field 
theory. They also proposed that an initially magnetized 
Fermi gas will tend to stay spin-coherent for long times, 
even in the presence of magnetic field noise that is naively 



expected to cause strong dephasing, provided the system 
is close to the transition into a ferromagnetic state [9]. 

Using an optical lattice and engineering the band 
structure to get flat (dispersionless) bands is another in- 
teresting route to achieving ferromagnetism. Such 'flat- 
band ferromagnetism' [10] has the advantage that the 
ferromagnetic state occurs at weak repulsive interactions 
and can be theoretically analyzed in a reliable fashion; 
however, an existing theoretical proposal along these 
lines involves working with fermions in the p-band of a 
honeycomb optical lattice [11] which is a significant ex- 
perimental challenge. A recent work [12] has considered 
the possibility of ferromagnetism for strongly interacting 
fermions in optical lattices and studied, within a phe- 
nomenological Landau theory, the energetics of possible 
spin textures (such as hedgehog states, domain walls and 
skyrmions) which might arise in a trap. 

In this paper, motivated by earlier work and by on- 
going experimental efforts, we revisit the problem of fer- 
romagnetism in a harmonically trapped two-component 
Fermi gas. We begin by using a "local density approx- 
imation" (LDA), sometimes referred to in the literature 
as the Thomas- Fermi approximation, to describe this sys- 
tem. Within the LDA, we find that the mean field re- 
lease energy of the trapped gas (as well as the potential 
energy and the kinetic energy component of the release 
energy) provides a simple, albeit indirect, diagnostic of 
the ferromagnetic transition. We find that the formation 
of nonzero local magnetization in the trap causes a sup- 
pression of the atom loss rate via three-body collisions. 
This suppression competes with the growth of the loss 
rate as the interactions get stronger, leading to a peak 
in the atom loss rate at an interaction strength which is 
very close to, but slightly beyond, the ferromagnetic tran- 
sition point. We then show how one might incorporate 
magnetization gradient (or 'surface tension') terms in or- 
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der to go beyond the LDA. Our energy functional is akin 
to an earlier phenomenological Landau theory for ferro- 
magnetisni in an optical lattice [12] but explicitly keeps 
track of the spatial dependence of all the Landau the- 
ory coefficients which arise from density variations in the 
trap. Using this extended energy functional, we study the 
energetics of various spin textures in the case where the 
net magnetization, which does not relax in these quan- 
tum gases, is assumed to be zero. This corresponds to 
choosing the initial population to be the same for both 
hyperfine species of fermions. We show, in this case, 
that a hedgehog configuration of the magnetization has a 
lower energy than an 'in-out' phase-separated configura- 
tion with a domain wall. A similar phenomenon has been 
predicted for fermions in optical lattices [12]. Finally, we 
turn to the effect of anisotropic trapping frequencies in 
a harmonic trap. While such anisotropics can be incor- 
porated by a trivial rescaling of coordinates in the LDA, 
this is no longer true in the presence of surface tension 
which leads to a breakdown of the LDA. (A breakdown 
of the LDA has been observed [13] and theoretically ad- 
dressed [14-16] in the context of polarized superfluids in 
highly anisotropic traps.) We use our extended energy 
functional to numerically study how the hedgehog state 
distorts upon going from a spherically symmetric trap 
to an anisotropic cigar-shaped trap. We conclude with 
estimates which indicate that the magnetic dipole inter- 
action between atoms can be neglected for ^Li and ^°K. 



II. FERROMAGNETISM WITHIN THE LOCAL 
DENSITY APPROXIMATION 



The local density approximation (LDA) for a trapped 
Fermi gas corresponds to simply assuming that the above 
results apply locally in the presence of a trap potential 
T^(R). The ground state energy of this trapped Fermi 
gas is then obtained by minimizing the energy functional 



E[{p,{R)}] = Jd'R 



+ l-(R)^p,(R)-^M,p,(R) 



(4) 



where {pa-{R)} denotes the density profile of both spin 
species [p-|-(R), p|(R)]. Here we have introduced two La- 
grange multipliers p^ which act as chemical potentials 
for the two spin species and serve to impose the con- 
straints J (PRpaiR) = Na- The separate constraint on 
each spin component arises from the assumption that the 
two spin components correspond to the lowest two Zee- 
man split hyperfine levels of a Fermi gas. Since the Zee- 
man splitting near a Feshbach resonance is typically far 
greater than all other energy scales and the total energy 
must be conserved in these thermally isolated gases, we 
arrive at the constraint that the population of the two 
Zeeman components cannot change for fermionic atoms 
where the only interaction is between different spin com- 
ponents. Thus, unlike in solid state ferromagnets, the 
magnetization can be conserved on very long time scales. 



The Hamiltonian describing a uniform two-component 
Fermi gas interacting through a repulsive s-wave contact 
interaction is given by 



H 



rCKC 



.+9 



Rcp^i c^j^ Cj^j^ Cj^i , 



where ek = ?i^K^/2M is the kinetic energy of atoms 
with mass M and momentum ?iK. For a Fermi gas with 
N„ particles of spin-o-, the uniform gas densities of each 
spin is Per = Na/V, and the total kinetic energy of the 
uniform gas is just 



K = -V{Ef^P^ + Ef^P^), 



(2) 
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where V denotes the system volume, and Epc^ = apj , 
with a = (67r2)2/3;i2^2M, denotes the Fermi energy of 
particles with spin-cr. A mean field theory of the inter- 
acting Hamiltonian yields the total interaction energy 



U = gVp^Pl. 



(3) 



At this level of treatment the contact interaction strength 
g is related to the two-body scattering length as in vac- 
uum via g = Airash^/M. 



Rescaling to the isotropic problem for a 
harmonic trap 



Let us assume that the trapping potential is harmonic, 
(1) but possibly anisotropic, so that ^(R) = iM^.w^^R?. 
Here 0Jx,y,z are the trapping frequencies along different 
directions. In order to make the problem appear isotropic 
we can rescale distances by setting Rj = Rj(a;j/r2), where 
O = {ijJxOJyOJzY^^ is the geometric mean of the trap fre- 
quencies. With this standard rescaling, we get 



E[{p,{R)}] = SR 



+ iMfi2R2^p^(R)_^^^p^(R) 



B. Noninteracting unmagnetized gas 



.(5) 



For the unmagnetized gas, we have TVt 



N/2 



and for the noninteracting case we can set g = Q. This 
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reduces the energy functional to 



(6) 



Setting dE*^ / dpa = leads to the equations 



(7) 



(8) 



where we have used symmetry to set = Atj, = The 
solution to this equation is simply 



1 



p.(R)=a-3/2 
Clearly there is a maximum radius 



3/2 



R% 



N 



(9) 



(10) 



beyond which p<t(R) = 0. Integrating upto this maxi- 
mum radius using spherical symmetry of the density, and 
employing the constraint, wc end up with 



E% 



= E% = hn{3N) 
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|f=aHo(24^)V^ 



47V 



(37V)4/3, 



(11) 

(12) 
(13) 



%S(^)'/'^'/^ (14) 



where auo = {h/Miiy^^ is the oscillator length, and 
we denote the density solution for this noninteracting 
unmagnctizcd Fermi gas by p^^(R). Here is the 
chemical potential of the gas, E*^ is the total energy of 
the gas, and p%{0) is the atom density at the trap center. 



C. Converting the interacting problem to 
dimensionless variables 



Let us use the noninteracting unmagnetized gas results 
to convert to dimensionless variables as follows. 



r 


R 

~ "pO"' 
-ttjV 


(15) 


Ha 


p%A^=oy 


(16) 


A 


= 4(0)as, 


(17) 


£ 


= E/E%, 


(18) 


7<7 




(19) 



Here A is the dimensionless interaction parameter, e, 7^ 
are the total energy and chemical potential respectively 
in dimensionless units, and kp{0) denotes the Fermi 
wavevector at the trap center for the unmagnetized non- 
interacting gas. 

In terms of these dimensionless variables, the energy 
functional becomes 



(nf (r)+nj'^(r)) 



5/3/ 



16 
3^ 

+^'^T(r)«i(r)-^(7<T-r^)na(r) 



• (20) 



If we assume that the ground state solution for the densi- 
ties respects the spherical symmetry of this energy func- 
tional, we can further simplify the energy functional to a 
one-dimensional integral 



^(nf (r)+n°/^(r)) 



5/3/ 



4A 



+^ {r)ni (r)-^ {la-r^ha (r) 



(21) 



D. Variational minimization 

The variational minimization 5E / 5n„(r) leads to the 
following two equations 



n|(r) = 

ni{r) = 



(7T -r"^ - -^^^^^'^)) 



3/2 



3/2 



subject to the constraints 
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47r / dr r^na{r) = 



4 N ' 



(22) 
(23) 

(24) 



These equations can be iteratively solved (numerically) 

for the fermion densities given the interaction strength 
and the total fermion numbers for each species. Having 
solved them we can use the resulting fermion densities to 
compute physical observables. We will denote the aver- 
age magnetization by rn = {Ni^ — Ni)/{N-^ + Ni). For 
an unmagnetized gas, we find that increasing the inter- 
action A progressively modifies the density profile of the 
gas from that of a noninteracting Thomas-Fermi profile 
at A = to that of a fully polarized gas when A 1. 



E. Release energy and 'ferromagnetic transition' 

The release energy of the trapped atomic gas is mea- 
sured by rapidly switching off the trap potential and mea- 
suring the total kinetic energy of the atoms after some 
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FIG. 1: (a): Dimensionless mean field release energy, £rci, 
versus interaction, A, within the LDA for indicated values of 
m = (A'^i — Ni)/N. For rh — 0, there is a phase separation 
transition which appears as a sharp kink in the release energy 
at A « 1.84. (b,c,d): Same as in (a) but for the kinetic energy, 
interaction energy and potential energy of the gas. The energy 
per particle in physical units may be obtained by multiplying 
these results by SEp/A where Ep = hQ{3Ny'^ . As shown 
later, going beyond the LDA leads to negligible quantitative 
corrections to these results. 



time delay. Assuming that the switch-ofF process is in- 
stantaneous and that all the interaction energy in the 
initial state has been converted into the kinetic energy 
of atoms at the time of measurement, the release energy 
and its separate kinetic and interaction energy compo- 
nents are given, within the LDA, by 



^roi = ^ Jdr r2 ^^n5/3(r)-)-i-Ant(r)n|(r 



— An|(r)n|(r) 
on 



£kii 



5^' 



,5/3 



(r) 



,(25) 
(26) 
(27) 



The potential energy of the cloud, due to the confining 
harmonic trap, can be easily obtained from experimental 
measurements of the cloud profile and it is given by 

Spot = ^Jdr r'[r'in^{r) + n^ir))] . (28) 

From Fig. [1] we sec that the release energy displays a 
sharp transition point, for m = 0, at Acrit ~ 1.84. At this 
interaction strength, we find that kpas = 7''/2 at the trap 
center, with kp being the Fermi wave vector at the trap 
center in the interacting cloud, which corresponds to the 
Stoner transition point in the uniform gas. Further, an 
examination of the density profile of the two spin species 
shows that this corresponds to an onset of phase separa- 
tion in the trap — for A > Acrit, atoms of one spin type 
tend to have a higher density at the trap center while 
atoms of the other spin type are pushed away from the 
center leading to a nonzero magnetization density near 
the trap center. Exactly which atoms tends to accumu- 
late at the center is a spontaneously broken symmetry at 
zero magnetization, and this phase separation is simply 
a local manifestation of ferromagnetic ordering. This re- 
sult for the Acrit at m = translates into an estimate for 
the critical two-body scattering length. 



(crit) 



O.eAcritOHO^ 



-1/6 



moN' 



-1/6 



(29) 



beyond which one expects to see phase separation in 
the trap. For n/2Tr « 170Hz we estimate for N = 
10^, 10^, 10^, the respective critical scattering lengths 



(crit) /40 



''S 

(crit) /g 



(^'^K) « 5300ao, 3600ao, 2500ao, (30) 
(^Li) w 13800ao, 9500ao, 6400ao, (31) 



where qq w 0.529A is the Bohr radius. 

Fig. [1] also shows that the kinetic energy and the in- 
teraction energy components of the total release energy. 
Each of these observables shows a large and much more 
dramatic signature at the transition (for m = 0) than the 
total release energy. It would be promising to look for 
this signature in experiments. In addition, the potential 
energy of the confined cloud also shows a maximum at 
the ferromagnetic transition point. 



5 



Strictly speaking, there is no phase transition (beyond 
mean field theory) except in the thermodynamic limit 
which, for a trapped Fermi gas, is obtained by taking 
N-^oo and fl-^O with Nil^ held fixed. For nonzero mag- 
netization, however, there is no phase transition even at 
mean field level; nevertheless, the release energy does dis- 
play a fairly sharp crossover at Acrit for m ^ 0.2. The 
measured release energy can only tell us about the ex- 
istence of a phase transition — for A > Adit, in situ 
measurements of the magnetization profile, which we dis- 
cuss below, are needed to show that this transition cor- 
responds to ferromagnetism in the trap. 
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F. Atom loss rate 

Atoms on the repulsive side of the Feshbach reso- 
nance tend to be unstable to formation of molecules via 
three-body collisions. Apart from kinematic and statis- 
tical contraints on these processes, there is a simple con- 
straint that two of these atoms, which eventually form 
the molecule, must have opposite spins. One consequence 
of having a nonzero local magnetization in the trapped 
gas is a suppression of the probability of finding fermions 
with opposite spin in the same region of space, which 
leads to a strong suppression of such three-body losses. 
A measured drop in the atom loss rate as a function of in- 
creasing interaction strength would thus hint at the pres- 
ence of nonzero local magnetization in the trap. Upto an 
unknown prefactor, Fq, we can estimate this three-body 
loss rate as 

r^FoA^y d^r nT(r)nx(r)(nt(r)+n^(r)), (32) 

where the A'* scaling follows from a study of the three- 
fermion problem [6]. Fig. [2] depicts a plot of F/Fq as a 
function of the interaction strength A. For to = 0, the 
very rapid growth of F/Fq for small interaction strength 
arises from the rapidly growth of the A^ coefficient, while 
the rapid drop beyond the ferromagnetic transition point 
arises from the formation of a nonzero magnetization 
which suppresses the product nf(r)ri|(r) in the inte- 
grand. These two competing effects lead to a peak in 
the rate of atom loss, via conversion to molecules, at an 
interaction strength which is slightly beyond the ferro- 
magnetic transition point. 



FIG. 2: Dimensionless atom loss rate, F/Fo, as a function of 
interaction strength at various average magnetizations. 



LDA by considering magnetization gradient terms with 
a view to studying the energetics of spin textures. We 
begin by noting that although we have been assigning a 
global spin axis to the magnetization, the LDA energy 
functional would be unchanged if we in fact choose the 
local spin quantization axis to vary from point to point; 
only the magnitude of the local magnetization plays a 
role. In order to go beyond the LDA and to study the 
energies of various spin textures in such a Fermi gas, we 
therefore need to extend the energy functional in two re- 
spects. First, we have to promote the local magnetization 
to a vector quantity so the magnetization can point in dif- 
ferent directions on the Bloch sphere at different spatial 
locations. Second, we have to include terms in the en- 
ergy functional which depend on the local magnetization 
gradients; this corresponds to adding a 'surface tension' 
term to the energy functional. The results from such an 
extended energy functional should be compared, in the 
future, with microscopic Hartree-Fock calculations. 

We start with the dimensionless energy functional in 
Eq.® and set 

fl( Y] 

^T(r) = -^(l + TO(r)), (33) 
n,(r) = !^(l-TO(r)), (34) 



III. BEYOND THE LDA: MAGNETIZATION 
GRADIENTS 

The discussion in the preceding section has focused on 
the properties of the Fermi gas within the LDA. The en- 
ergy functional at this level of approximation does not 
have any gradient terms. We will not worry about the 
shortcomings of this approximation for the density pro- 
file — it is well known that the LDA breaks down near 
the trap edges — but instead focus on going beyond the 



which defines the local magnetization density m(r). As 
discussed, the spin quantization axis can be chosen to be 
different at each space point within the LDA. Let us next 
expand the energy functional in powers of TO(r); we will 
keep terms upto m^(r) although terminating the expan- 
sion at m^(r) would not qualitatively affect our results. 
The energy functional then splits into two parts as 

e = ea[n{r)]+ et[n{r),m{r)], (35) 



6 



where 



£a[n(r)] = ^ 



6."(r)N5/3 
5^ 2 



eb[7i(r),m(r)] = 



16 

(7 - r2)n(r)] , 

i^|rf^r[A,(r)m^(r) 

A4(r)m4(r)+A6(r)m6(r) 
ft, n(r)m(r)] . 



(36) 



(37) 

Here only depends on the density profile which de- 
pends on the interaction A and which we assume is un- 
changed from that given by the LDA calculation earlier. 
This is a good approximation since the corrections to the 
LDA energy are weak for typical atom numbers used in 
experiments as we will see below. The coefficients of the 
magnetization-dependent energy functional, et, are 



Mr) 



^ 22/33 
n5/3(r) 
22/381 ' 
7n5/3(r) 



^-^(r)), 



(38) 
(39) 
(40) 



22/32187' 

A2{r) depends on A explicitly. In addition, all the coef- 
ficients ^2,4,6 depend on the spatial location in the trap 
through the density, and thus also depend implicitly on 
the interaction strength A. This dependence was ignored 
in earlier phenomenological work on trapped fermions in 
an optical lattice [12]. The Lagrange multiplier in the 
energy functionals are given by 7 = (7^ -I- 7x)/2 and 
h = (7f — 7x)/2. Promoting the magnetization and the 
Lagrange multiplier h to vectors m, h, and including gra- 
dient terms leads to an energy functional of the form to 

/d3r[A2(r)|m(r)|2 



eb[n{r),m{r)] 



37r2 



-A4(r)|m(r)|4 + A6(r)|m(r)|« 
-iCs(i")aj(Vimj(r))^-K(r)-m(r) 



,(41) 



Here, Ui = (uji/il)'^ which comes from our rescaling to 
an isotropic problem. The stiffness Cs(r) depends on r 
only through the density n{r), and it can be computed in 
the uniform Fermi gas assuming that the magnetization 
variation is slow on the scale of the intcrparticle spacing, 
but fast on the length scale over which the total density 
varies, so that density variations can be ignored in this 
computation. The Lagrange multiplier h{r) must be cho- 
sen to satisfy global constraints on the magnetization, for 
instance, / d'^r n(r)mi(r) — for each component i. We 
next outline the derivation of the stiffness term. 



A. Computation of the stiffness Cs(r) 

For small magnetization, we can obtain the result 
for from the result for the magnetic susceptibility 



of the uniform Fermi gas. Note that the excess en- 
ergy in an applied field ft,(q) (pointing in any direction) 
is given by A£'(q) = ^x(q)^i(q)^j(— q) which defines 
the wavevector dependent magnetic susceptibility. This 
tells us that the magnetization M(q) in this external 
field is simply M(q) = x{q)h{q), so that we can in- 
stead set A£;(q) = ^x~^ {<l)Mi{q)Mi{-q) . Expanding 
Xq^{^ + bq^) then yields 



X-^(q 



Ai?(q) - -xo 1(1 + 6q2)M,(q)M,(-q) 



(42) 



The well-known result for a Fermi gas at T = is that 
b = l/12kp, using which the energy cost becomes, in real 
space, 



AE-- 



where 



2X0 



1 



|A^(R)r + Y2fcr(VA/,(R)) 



Xo 



Mkp 



■(3^V) 



^1/3 



(43) 



(44) 



Rescaling distances to get an isotropic harmonic trapping 
potential, and setting Mi = p%„{r — 0)n{r)mi{r), with 
r = R/R%, we find 



Cs{r) 



^1/3, 



r) 



1 



22/33 6(37r2n(r))2/3 \p%^{R%)^ 
1 



1 



2/3 



72n(r)(3iV)2/3- 



(45) 



For general values of the magnetization, higher order gra- 
dient terms might also become important. We will focus 
here on the effects of this simplest gradient term in the 
energy functional. 



B. Simplified magnetization energy functional 

Before proceeding to the energetics of various spin tex- 
tures, let us slightly simplify the energy functional. No- 
tice that n(r) varies over the length scale of 1 in our di- 
mensionless units. For large atom numbers, the stiffness 
is small as seen from Eq. psj) and we therefore expect 
significant variations of the magnetization to occur on 
length scales £ ^ 1 in our dimensionless units. Making 
this assumption, we can set V(rt(r)mi(r)) « n(r)Vmi(r), 
which results in the slightly simplified energy functional 



£6 



^ J d'r [A,{r)\m{r)\' + A,{r)\m{rt 



+ A6{r)\m{r)f + h{r) ■ m{r) 
n{r) 



+ 



144(3iV)2/3 



where h{r) is chosen to satisfy 



y d^r n{r)mi{r) — 0, 



(46) 



(47) 
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for each component i (for zero net magnetization). Re- 
call that ai = [uji/^Y, where = (uj^UyUJzY/^ is the 
geometric mean of the trap frequencies. 



IV. ENERGETICS OF SPIN TEXTURES 

The energy functional we have derived above allows 
us to study the energetics of various magnetization pat- 
terns in the trapped Fermi gas. We begin by consider- 
ing the isotropic harmonic trap, for which we compare 
energies of a hedgehog configuration and a domain wall 
configuration of the magnetization. We then turn to an 
anisotropic cigar-shaped trap and show how the hedge- 
hog state gets deformed from the isotropic case. In each 
case, we begin by constructing the appropriate ansatz for 
the magnetization. We then numerically minimize the re- 
sulting energy functional, by discretizing it on a fine grid 
of points and using a simulated annealing procedure, to 
obtain the optimal magnetization profile and its energy. 
Assuming that the density and magnetization satisfy the 
constraints that the total atom number is fixed and the 
total magnetization is zero (so that the Lagrange multi- 
pliers can be dropped), we can express the total energy 
as a sum e = ei + £2 where 



£1 



16 
3^ 
16 



r n(r) 



, d\{A2{v)\rh(v)f ^ A^(v)\Tn(j)\' 



+ A(r)|m(r)|6 + 



rt(r) 



144(3iV)2/3 



ai(V,mj(r))^ 



(48) 



.(49) 



A. Isotropic trap: Hedgehog state 

For the isotropic trap, the density profile is spherically 
symmetric, which allows us to set 



64 2 
37r 



5 2 6t\ 



(50) 



For the magnetization-dependent energy functional, we 
must set ai = 1 in the isotropic trap, and the hedgehog 
state corresponds to choosing m(r) = m{r)f. This leads 
to 



£2 



64 
3^ 



dr [j42 (r)m^(r) 



-I- Ai{r)'m'^{r) + Afi{r)m^{r) 
n{r) 



144(37V)2/3 I r-2 



dm{r) 
dr 



.(51) 



We do not have to pay attention to the zero magnetiza- 
tion constraint since this is guaranteed for any choice of 
m{r) by the hedgehog ansatz symmetry. For typical par- 
ticle numbers in experiments, N ^ 10^ — 10^, the stiffness 
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FIG. 3: (Color online) Dimensionless total energy, ei +82, 
shown as a function of interaction strength. A, for an isotropic 
harmonic trap. DW indicates the energy of the domain wall 
state for 10* atoms (solid) and 10® atoms (dashed). HH de- 
notes £1 -I- £2 for the hedgehog state which is nearly identical 
for 10* and 10® atoms. Also shown (thin solid line, 'unmag- 
netized') is ei, defined in Eq.(44), which depends only on the 
total density profile. 



term has a very small coefficient. We will therefore as- 
sume that the average density profile n(r) obtained from 
our earlier LDA calculation remains unchanged and only 
focus on changes in the magnetization profile arising from 
inclusion of gradient terms. 

Fig. [3] shows the energy e^^ of the hedgehog state ob- 
tained by finding the function m(r) which minimizes the 
hedgehog state energy. Fig.[Ua) shows the magnetization 
profile of the hedgehog state at two different interaction 
strengths. We find that the magnetization is suppressed 
in a small region around the trap center and vanishes at 
r = 0. To understand the magnetization profile of the 
hedgehog near the trap center, we can focus just on the 
last two terms in Eq. ([5T|) . Taking a functional derivative 
with respect to m{r) and setting it to zero then suggests 
that m(r) ^ at small r, so the energy density com- 
ing from the central region of the hedgehog is finite. Far 
from the center, we expect the magnetization to be small. 
These expectations are consistent with the magnetization 
profiles shown in Fig. HJa). 



B. Isotropic trap: Domain wall state 

For the domain wall state we have, as before, 

Xn\r) 



64 / , 2 
El = — fdr r 

on 



6 ,n{r) 
5^^ 



(52) 



For the magnetization dependent energy functional, we 
set ai = 1 in the isotropic trap and choose m(r) — m(r)z 
for the domain wall. This is capable of describing a state 
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FIG. 4: (a) Magnetization profiles for the hedgehog state at 
A = 2.0 (solid) and A — 2.4 (dashed), (b) Magnetisation pro- 
files for the domain wall state at A = 2.0 (solid) and A — 2.4 
(dashed). The profiles have been calculated for 10^ atoms in 
an isotropic trap. The hedgehog state has zero magnetization 
at the trap center while the domain wall state magnetization 
gets suppressed around the domain wall but remains nonzero 
at the trap center. Insets indicate the schematic magnetiza- 
tion plot of the hedgehog state and the domain wall state. 



with spin-t, say, at the trap center with spin- J, pushed 
away from the center, what we might call an 'in-out' 
domain wall. We find 



£2 



64 
3^ 



dr [j42(r)TO^(r) + A4^{r)'m'^{r) 



n{r) fdm{r) 



144(3iV)2/3 



dr 



(53) 



where, for iV| — N^, we must satisfy the constraint 
J dr r^n{r)m{r) = 0. Fig.[3]shows the energy e-^^ of the 
domain wall state obtained by finding the function m(r) 
which minimizes its energy subject to the zero magneti- 
zation constraint. We find that this domain wall state 
has a higher energy than the hedgehog state. Fig. |D^b) 
shows the magnetization profile of the domain wall state. 
As expected, the magnetization is suppressed in a small 
region around the domain wall but remains nonzero at 
the trap center. 



C. Cigar-shaped trap: Distorted hedgehog 

If we consider a cylindrically symmetric (cigar-shaped) 
trap, we can look at an ansatz of the form 

I X y \ 

m(r) — m{p,z) — sin ?/;, — sin ?/;, cos ip , (54) 
\P P J 



where p = ^ + y'^ and ip = ipip^z). For ip = 6 = 
tan^^(p/z) this reduces to the spherical hedgehog ansatz. 



Note that the direction of the magnetization on the Bloch 
sphere is unrelated to the location in real space. We could 
equally well have chosen, for instance. 



X y 
cos -0, — sin V', — sin -0 
P P 



m(r) ~ rn{p, z) 



With the choice of magnetization in Eq.([5 
|m(r)p = m'^ip, z), while 

{dim J ) {d^ruj ) = a± {dpmY + {d^mf 



(55) 
we have 



-— sin V 

' P^ 



rr? [ai(9p0)2 + a,(9,V)' 



(56) 



so the integral / S'v 27t J dzdpp. We can assume that 
TO is an even function of z and that ili{p, —z) ~ tt~'iIj{p, z) 
(so that sin^ 7/'(p, — z) = sin^ tp{p, z)) to restrict the en- 
ergy integration grid to just z > 0. These conditions 
ensure that the total magnetization integrates to zero. 
The final expression for the energy can thus be recast, 
with r = \/p2 + z^, as 



£2 

F{p,z: 



64 
3^ 



'dr r' 



6,?T.(r),r/Q Xn?(r) n , ,1 



64 



R„ 



dz 

A^rri^ 



dppF{p,z) (58) 



■ sin^ -0 



144(3iV)2/3 p- 
+ to2 [a^idp^py + a^id^^y]} , 



(59) 



with ip^p = 0, z) = and ^p{p,z = 0) = 7r/2 by sym- 
metry. For notational simplicity, we have suppressed the 
coordinate labels on n, to, tJj in the above functional. 

We find, numerically, that tp fa 9, so in fact the ansatz 
simplifies to the form 



m(r) — m(p, z) — , — , — 
\r r r I 



(60) 



The main effect of going from the spherical to the cigar 
shaped trap is that the magnitude of the magnetization 
is no longer just dependent on the radial coordinate r. 
The magnetization however still points (in our rescaled 
coordinates) along the radial direction. The plot of the 
magnetization for A = 2.4 in the rescaled and in the 
original coordinates for a trap anisotropy corresponding 
to aj_ = 2,q;z — 0.25 (a trap frequency ratio lo^/uJz ~ 
2.8) is given in Fig. [5l 



V. EFFECT OF DIPOLAR INTERACTIONS 

Our results for the spin texture energetics and mag- 
netization profiles have been obtained by neglecting the 
role of the long range magnetic dipole interaction be- 
tween the fermions. The dipole interaction will add to 
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FIG. 5: (Color online) Two-dimensional magnetisation pro- 
file for the distorted hedgehog showing breakdown of the LDA 
for the magnetization density for 10^ atoms in an anisotropic 
trap with A = 2.4, and a± =2,olz = 0.25 {uj±/uJz ~ 2.8). (a) 
Plot of the magnitude of the magnetization m(p, z) and equal- 
magnetization contours displayed in rescaled coordinates in 
which the trap potential is spherically symmetric. Colorbar 
to the right indicates the value of m(p, z). We see that m(p, z) 
is larger in magnitude for larger values of p than it is for z, in- 
dicating that the surface tension makes it easier to change its 
value in the weak trapping direction, (b) m shown as a quiver 
plot indicating the magnitude and direction of the magneti- 
zation (plotted in coordinates where the trap anisotropy is 
explicitly shown). Shaded area indicates the region of the 
trap where the atom density is nonzero. 



the magnetic energy of atoms in the trap. In addition, it 
will lead to spatial variations of the magnetic field seen 
by atoms within the trap and, thus, cause tend to cause 
dephasing as atoms in different regions will precess at 
different rates. Such effects are known to be important 
in *^Rb spinor Bose condensates [17,18]. In order to es- 
timate the dipole interaction energy and the timescale of 
this dephasing, we have considered the spatial variations 
of the dipole field for the simple case of the spherical trap. 

The expression for the precession frequency at distance 
r from the center of the spherical trap is, for the hedgehog 



dip / \ 1 Mo 2 / 



2%/2 3/2 



dridO r\ sin0 m{ri)n{ri) 
(r^ + — 2rri cos6')3/2 



i^(r,ri,0),(61) 



cos 6* — 3 



(ri — r cos 0) (ri cos — r) 



(r2 



2rri cos( 



(62) 



where jiB ~ 9.27x is the Bohr magneton, and 

/io = 47r X 10~^7V/A^ is the permeability of free space. 
Evaluating this, we find that the typical value of (and also 
the variation in) the precession frequency for A — 2.5, 
N = 10^ and VL = 27r(170Hz), is i^^']^(^Li) w 0.03Hz 
and v'^'^C^^K) « 0.6Hz. The energy associated with the 
dipole interactions is far smaller than our estimate of 
magnetic exchange energies, ~ SOOHz, arising from the 
s-wave contact interaction between fermions (in the in- 
teraction range where we expect ferromagnetism) . At the 
same time, measurements of the typical atom lifetime, r, 
on the repulsive side of the Fcshbach resonance indicate 
that T - 10ms for '^^K [19] and r - 100ms for ^Li [2,20]. 
These are clearly much less than the variations in the pre- 
cession period induced by spatial variations of the dipolar 
field as estimated above. Taken together, these estimates 
show that ignoring the effect of dipole interactions is a 
very good approximation in this system. 



VI. CONCLUSIONS 

In conclusion, we have studied ferromagnetism and 
spin textures in ultracold atomic Fermi gases in the 
regime of strongly repulsive interactions using the LDA 
extended to include magnetization gradient corrections. 
Within the LDA at zero temperature, we have shown 
that the release energy of the gas, as well as its sepa- 
rate kinetic energy and interaction energy components, 
shows a sharp signature of the ferromagnetic transition. 
We have also shown that the atom loss rate via three- 
body collisions has a peak very close to the ferromagnetic 
transition and it provides yet another diagnostic of the 
transition into the ferromagnetic state. We have gone 
beyond the LDA by deriving a surface tension correction 
to the energy functional, which depends on atom number 
and the trap-geometry, and used it to study the energet- 
ics of various spin textures in a two-component trapped 
Fermi gas. For a spherically symmetric trap, we find 
that a hedgehog magnetization profile has lower energy 
than a domain-wall state. For large atom numbers, the 
small surface tension leads to a small energy difference 
between the two spin textures and the results are close 
to those of the LDA. In this case, the surface tension is 
responsible for selecting the hedgehog state as having the 
lowest energy but we have checked that it does not sig- 
nificantly change our results for the release energy and 
the atom loss rates. These continue to be useful, albeit 
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indirect, diagnostics of the transition into the ferromag- 
netic state. For elongated clouds, we have shown that 
the surface tension term distorts the hedgehog states, in 
rescaled coordinates where the trap is isotropic, in such a 
manner that the magnitude of the magnetization changes 
more easily in the weak direction of the trap than would 
be expected on the basis of the LDA. Such a breakdown 
of the LDA is more apparent for smaller atom numbers. 
Finally, we have considered the effect of magnetic dipolar 
interactions on our results and find that it is a good ap- 
proximation to ignore dipole interactions in this system. 

The typical atom loss rate near the Feshbach resonance 
sets a constraint that the formation time for the ferro- 
magnetic state will have to be on the order of tens of mil- 
liseconds for ''"K [19], and hundreds of milliseconds for 
^Li [2,20], in order for it to be observed. A direct way to 
probe the spin textures discussed here would be through 
high resolution in situ magnetometry as has been done 
for spinor Bose condensates [21]. An experimental obser- 
vation of ferromagnetism in trapped Fermi gases would 
provide impetus for future theoretical work on finite tem- 
perature effects and collective modes in the strongly in- 
teracting regime. 
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